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Abstract 

Considering a basic enzyme-catalysed reaction, in which the rate of 
input of the substrate varies periodically in time, we give a necessary 
and sufficient condition for the existence of a periodic solution of the 
reaction equations. The proof employs the Leray-Schauder degree, 
applied to an appropriately constructed homotopy. 



1 Introduction 



The basic scheme for a reaction catalysed by an enzyme is 

I(t) h k 2 

-> S + E +± C -> P + E 
fc-i 

in which the substrate S and the enzyme E form a complex C through a 
reversible reaction, and the complex C can dissociate into the enzyme and 
the product P. I(t) is the rate of input of the substrate into the system, 
satisfying 

I{t) > 0, tel. (1) 
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The dynamics of this system are described by the rate equations for the 
concentrations of the species: 



E\t) = -hE(t)S(t) + (A;_i + k 2 )C(t) (2) 

S'(t) = I(t) - k x E(t)S{t) + fc_ x C(t) (3) 

C'{t) = hE(t)S(t) - (A;_i + A; 2 )C(t) (4) 

P\t) = k 2 C(t) (5) 



Stoleriu, Davidson and Liu [2] recently investigated the case in which the 
rate of input of the substrate, I(t), fluctuates in time in a periodic manner. 
As they noted, such a situation is common in biological systems, both due 
to intrinsic oscillations in preceding steps of the reaction pathway, and to 
oscillations external to the organism. 

Since by adding © and (TJ) we have [E(t) + C{t)}' = 0, that is the total 
amount of the free and bound enzyme is constant in time, we can set 

E(t) + C(t) = K, 

and since (jSJ) decouples from the other equations, we can rewrite the system 
in terms of S and C as 

S'(t) = I(t) -h[K- C(t))S(t) + fc_iC(t), (6) 

C'{t) — ki[K — C(t)]S(t) - (fc_! + k 2 )C(t). (7) 
Assuming that I(t) satisfies §I§ and is T-periodic 

I(t + T) = I(t), teR, (8) 

we ask, as a first step in understanding the dynamics of the system, whether 
there exists a T-periodic solution S(t), C(t), of (jBJ), ©, with 

S{t) > 0, ieR, (9) 

< C(t) < K, teR. (10) 

In the case of a constant rate of input I(t) = I, there is a unique stationary 
solution if and only if 

I<k 2 K, (11) 
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given by 

c ~k 2 h s - kl (k 2 K-iy (12) 

and a phase-plane analysis shows that all solutions tend to this equilibrium. 
Note that the condition ffTTj) says that the rate of input I of the substrate 
is not too large - indeed if this condition is violated then S(t) will increase 
without bound, since the substrate enters the system more rapidly than it 
can be processed by the enzyme. The proof of our existence theorem for the 
periodic case will involve a homotopy connecting it to the case of constant 
input. 

Returning to the general case of periodic I(t), and adding and ([7]) we get 

C'(t) + k 2 C(t) = I(t) - S\t). (13) 

Assuming that C(t),S(t) is a T-periodic solution, and integrating ([TBI on 
[0,T], taking into account the periodicity, we get 

k 2 [ C{t)dt = [ I(t)dt. (14) 



From ([IU]) and (dD it follows that 

I(t)dt<k 2 K. (15) 

Thus fTTS]) is a necessary condition for the existence of a T-periodic solution. 

In [2J, (where the case I(t) = 7(1 + esin(o;t)) with < e < 1 is considered) 
it is proposed that the existence of a periodic solution can be proven by 
constructing an invariant rectangle for the flow corresponding to the system 
©, ©, of the form 

D = {(S,C) \0<S<S, 0<C<K-5}, (16) 

(with S and 5 chosen appropriately), so that Brouwer's fixed-point theorem, 
applied to the time T Poincare map of ([51), (171), implies the existence of a fixed 
point of this map, which corresponds to the required T-periodic solution. An 
examination of this method of proof shows that an invariant rectangle of the 
form (1161) exists if and only if the input function I(t) satisfies 

max lit) < k 2 K. (17) 

te[o,T] 
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Note that this sufficient condition (IT7|) for the existence of a T-periodic solu- 
tion is more stringent than the necessary condition (j!5p . Here we bridge this 
gap by proving that in fact ffl5|) is sufficient for the existence of a positive T- 
periodic solution. We employ the methods of nonlinear functional analysis, 
reformulating the problem as a fixed-point problem for a nonlinear operator 
in a space of T-periodic functions, and using degree theory (see, e.g., [HE]) 
to prove the existence of a fixed point of this operator. 

Our existence result raises the following question: is it true that for any I(t) 
satisfying ffl5l) . the periodic solution is unique and globally stable! 

2 The existence theorem 

We set 

I='~J I(t)dt, I (t) = I(t)-I. 

Theorem 1 Assume k±,k2 > 0, k-\ > 0, K > and I(t) is a continu- 
ous function satisfying (TJP and (TJJ]. Then there exists a T-periodic solution 
S{t),C{t) of satisfying if and only if 

< I < k 2 K. (18) 

The fact that (1181) is a necessary condition for existence was explained above, 
so we now assume that ffTB"]) holds, and need to prove the existence of a T- 
periodic solution. 

Note that © , © is equivalent to (|6|) . (fl3l) . and it is the latter which we will 
be using. 

We define 

0=1-1, (19) 

Q(t) = C(t) - 0. 
The condition fflBl is equivalent to 

< C < K, (20) 
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and from (]T4"j) we have that if S, C is a T-periodic solution of (l6"j). (lT3"]) then 
C satisfies 

r-T 

C (t)dt = 0. (21) 



We can thus rewrite (161. (113]) in terms of S and Co: 

S'(t) + h[K - C - C (t)]S(t) = I(t) + fe_iC (t) + fc-iC, (22) 

C (t) + A: 2 Co(t) = /o(t)-5 / (t). (23) 

The following lemma, in which we solve (1221) for a T-periodic solution S(t) 
in terms of C (i), is the place where the key condition (1TB]) is exploited. 

Lemma 1 Given a T-periodic continuous function Co(t) satisfying l[21\) . the 
linear equation [2^1 has a unique T-periodic solution S(t), and if, in addition, 

C(t) = C + C (t) > 0, iel, (24) 

then S satisfies (Tj|). 

Proof: Setting C{t) = C + C (t), the general solution of ([22]) is 

S(t) = exp (- fex y [K -C(s)]ds^S(0) (25) 
+ y exp (An J [K- G{r)]drj [I(s) + k^C{s)]ds. 
The condition for 5 to be T-periodic is S(0) = S(T), or 

1 - exp^-^y - C(s)]ds) 5(0) (26) 
exp (hi y [X - C(r)]dr) [I(s) + fc^C^t/a. 



^iT^-C] 



(120)) , (ITS]) and (1211) imply that 



f [K- C{t)]ds = T[K - C] > 0, 
Jo 



which implies that the coefficient of S(0) on the left-hand side of (l2T)j) is 
positive, hence the existence of a unique T-periodic solution S of (f2"3j) is 
ensured, given by 

S(t) = J exp (ki J [K~ C(r)]dr) [I(s) + fc_iC(s)]ds (27) 

+ [^k-G) - l]- 1 J* exp (kt J\k -C{r)}dryi{s)+k^C{s)]ds. 

If (]2l]) holds, then recalling (JTj) and (1201) . we see that the second term in (1271) 
is strictly positive, and since the first term is nonnegative, we get (J9]). ■ 

We define X to be the space of continuous T-periodic functions C (t) satis- 
fying fl2"Tj) . with the maximum norm. 

Fixing I(t) satisfying ATS]) , we use the result of lemma [T] to define a mapping 
^(J; ■) : X — ► X (J is regarded as a parameter) as follows. Let Co G X 
and let S be the T-periodic solution of (1221 (given explicitly by (1271) ). Then 
define 

9 r (/;C ) = y. 

Note that since S is periodic the integral of 5" over [0, T] is 0, so we indeed 
have 9^(1; Co) G X. From the explicit formula (1271) it can be be shown by 
standard methods that 

Lemma 2 The mapping 5F(J; .) : X — > X is Frechet differentiable, and maps 
bounded sets to bounded sets. 

We can now reformulate the problem of finding periodic solutions of (|22|) . ( j23l) 
as: find C G X satisfying 

C' (t) + k 2 C (t) = I {t) - C )(t). (28) 

We define a linear mapping £ : X — > X: for any R & X, define = C 

to be the unique T-periodic solution of the equation 

C\t) + k 2 C{t) = R{t). (29) 

By integrating (|29~I) over [0, T] we see that 

/ C{t)dt = 0, 
6 



so that £j(R) G X. Moreover, since in fact L maps X boundedly to the space 
of C 1 -functions, which is compactly embedded in X, we have that 

Lemma 3 £ : X — > X is compact. 

We can now rewrite (I28I) as 

Co = L{I -V{I; Co)). (30) 

We will prove the existence of a solution of fl30l by applying degree theory. 
Fixing /, we define the homotopy !K : [0,1] x X — ► X by 

X(\,C ) = H(\I -3 ! (I + \I Q ;C )), 

and we will consider the equation 

Co = ^(A,C ) 0<A<1. (31) 

Note that for A = 1 (131]) coincides with (!30j) , while for A = (1311) corresponds 
to the system with constant rate I of substrate input. 

From lemmas [2] and [3] we have 

Lemma 4 : [0, 1] x X — > X is a continuous compact mapping. 

This ensures that we can apply the Leray-Schauder degree (see, e.g., [HE]) 
to IK. 

We define the bounded open set G C X by 

G = {C G X | < C„(t) + C < K t G E}. 

The following lemma summarizes the reduction of our problem to a fixed- 
point problem, and provides essential a-priori bounds. 

Lemma 5 If C G C solves / f57]) /or some < A < 1, and if C = C + C 
and S is the T -periodic solution of 

S'(t) = I + XI {t) -h[K- C(t)]S(t) + k^C{t), (32) 

C'{t) =k x [K- C(t)]S(t) - (fc_! + fc 2 )C(i). (33) 
and 5", C satisfy /fPj)./ f70j) . 
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Proof: Assume C satisfies ( T5TT) . By the definition of 3 we have 9^(7 + 
A/ ; Co) — S' where S is the periodic solution of (I3"2j) . and by the definition 
of £ we have 

C (£) + k 2 C (t) = XI (t) - 2(1 + AJ ; C ) = A/ (£) - S'(£). (34) 

Taking the difference of (I32I) and (1541 gives (|33|) . 
By ([T]) we have 

I + \I (t) = (1- X)I + XI(t) > 0, A £ [0,1], t El, (35) 
and the assumption that Co G G means that 

< C(t) <K, £ G M, (36) 
so lemma [1] implies that S satisfies Q). 

To show that C(t) > for all t, let t± be a point where C achieves its global 
minimum, and assume by way of contradiction that C(£i) < 0, so that by 
(I36p C(£i) = 0. Since t\ is a minimum point, we also have C'(t\) = 0. 
Substituting into (|3"3"|) we get kiKS(ti) = 0, so S'(ti) = 0, which contradicts 
(Q, which has already been proved. 

Similarly, to show that C(t) < K for all t, let £2 be a point where C achieves 
its global maximum, and assume by way of contradiction that C(£2) > K, 
so that by (1361) C(£a) = K. Since £2 is a maximum point, we also have 
C'(£ 2 ) = 0. Substituting into (13"3l gives + ^2)-^ = 0, a contradiction. ■ 

We note that lemma E] implies that if Co G C then Co satisfies ffTU]) . so that 
C G C. We thus get 

Lemma 6 A G [0, 1], C G <9G C ^ IK(A, C ). 

By the homotopy invariance property of the degree, it follows that 

Lemma 7 deg(id x — IK(1, •), C) = deg(id x — Ji(0, ■), G). 

We now show that 

Lemma 8 de#(id x - JC(0, ■), G) ± 0. 
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Proof: When A = 0, the only T-periodic solution of (1521) . (1551) is that given 
by (lT2j) . so the only solution of (15T1) is C = 0. 

In order to prove the lemma it suffices, then, to prove that the local degree 
of idx — IK(0, •) at Co = is nonzero, and this will follow if we can show the 
Frechet derivative idx — Dc 3t(0,Q) is nonsingular. 

We have 

id x - DcMl 0) = id x + £ o D Co 3(I- 0), 
so if this is singular there exists a nontrivial Co G X with 

C o = -L(Dc o ni;0)(Co)). (37) 

We will thus assume that (157|) holds, and show that it forces Co = 0. 
By the definition of the Frechet derivative, 

D Co ^(I;0)(C ) = ^(I; a C ) . (38) 
da a=o 

By the definition of 5F, we have 

V(I;aC ) = S t (a,t), (39) 

where S(a, t) satisfies 

S t (a, t) + h[K-C- aC (t)]S(a, t) — I + k-iaC (t) + k^C (40) 
In particular 

S(0,t) = S, 

where S is given in (1T2|) . Differentiating PUI) with respect to a and setting 
a = we get 

S Qt (0,t) + h[K- C]S a (0,t) = [kiS + k-t]C {t) (41) 
From (15H1) and (139|) we have 

D Co 3 r a;0)(c ) = ^(o,t), 

which, using (J57]l . implies 

£(5 Qt (0,-)) = -C , 
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so by the definition of £ 

C' (t) + k 2 C (t) = -S at {0,t) (42) 
Differentiating (|4ip with respect to t we have 

S aU {0,t) + h[K- C]S at {0,t) = [hS + k^]C' (t) 
and substituting (J4"2l) we have 

^«tt(0,t) - fci[# - + fc 2 (7 (t)] = + fc_J(75(t). (43) 

Differentiating (|42l) with respect to £ we have 

Satt(0,*) = -<%'(*) -*&<%(*), 
and together with (1431) we get 

Co (*) + [ fc 2 + h(K -C + S) + fc_i]C£(t) + £;iA; 2 [K - C]Q,(t) = 0. 

Multiplying this by Cg(t) and integrating over [0, T], taking into account the 
periodicity, we obtain 

[k 2 + h{K - C + S) + fc_x] / (C' (t)) 2 dt = 0, 

Jo 

and since, using (12"U|) . the coefficient is positive, we have C' = 0. Since 
C E X, this implies C = 0, as we wanted to prove. ■ 

From lemmas [7] and [S] we obtain 

deg(id x -K{l,-),G) ^ 

which implies the existence of a solution Cq E G of (l3Tj) . which, by lemma [5j 
implies theorem [TJ 
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